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• We study a one-dimensional fluid of hard-rods interacting each other via binary inelastic collisions 

' and a short ranged square-well potential. Upon tuning the depth and the sign of the well, we 

investigate the interplay between dissipation and cohesive or repulsive forces. Molecular dynamics 
simulations of the cooling regime indicate that the presence of this simple interparticle interaction 
l/") ' is sufficient to significantly modify the energy dissipation rates expected by the Haff's law for 

the free cooling. The simplicity of the model makes it amenable to an analytical approach based 
on the Boltzmann-Enskog transport equation which allows deriving the behaviour of the granular 
temperature. Furthermore, in the elastic limit, the model can be solved exactly to provide a full 
thermodynamic description. A meaningful theoretical approximation explaining the properties of the 
. inelastic system in interaction with a thermal bath can be directly extrapolated from the properties of 

the corresponding elastic system, upon a proper re-definition of the relevant observables. Simulation 
results both in the cooling and driven regime can be fairly interpreted according to our theoretical 
approach and compare rather well to our predictions. 
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' PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 
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I. INTRODUCTION 



e 

Granular materials are ubiquitous in nature and their handling occurs in many types of industrial activities. While 
they are very common, their properties often are not. In the last twenty years there has been a great progress in 
the comprehension of static and dynamical properties of granular fiowsi&Si^. In spite of the fact that most of the 
theoretical research in this context has been based on the inelastic hard sphere model, several observations suggest that 
I ' neither cohesive forces^ nor electrostatic repulsionZ can be ignored. Understanding how simple interactions modify 
. the behavior of a granular gas can have important practical consequences. Cohesive forces have to be considered 
when studying wet granular matter: the humidity may lead to the formation of thin layers of water on the surface 
of the grains and induce adhesion through capillarity effects. The presence of liquid-vapor interfaces can enhance 
1 the mechanical stability of an assembly of grains, as illustrated by sand castles 8 . On the other side, repulsive forces 
also play a role, as stressed by Sheffler and Wolf 7 . Dry granular materials tend to become electrically charged due to 
■ contact electrification during transport. In the case of monopolar charging the particles experience mutual coulombic 
repulsion. Finally, Blair and Kudrolli 9 studied the behavior of a vibrated system of magnetic grains, where forces 
of tensorial character are in action, and found coexistence of long lived clusters with isolated particles. Clusters can 
' manifest as chains or globular structures according to the driving intensity. 

In this work we introduce and study a one-dimensional model which can be tuned to describe both the cohesive and 
the repulsive regime. One dimensional models have often been employed in the literaturoi2iii*i2i±iii4 to study granular 
fluids because their simplicity provides a valuable testing ground for theoretical approaches and approximations. Our 
model consists of a set of inelastic hard rods subjected to square-well potential as shown in Fig. The attractive 
potential mimics the action of cohesive forces responsible for adhesion among particles which are crucial effects when 
considering fine particulates such as powders or sands. On the contrary, the barrier describes the effect of soft materials 
which may present a deformable shell covering the hard-core nucleus. 

The choice of a square well interparticle interaction is particularly convenient in a computer implementation of 
the model since it reduces Newton's equations to algebraic expressions. Indeed, in the cooling regime, rods move 
with constant velocity until they pass a barrier or their cores touch. Thus the collisional cooling can be simulated 
through the collision driven algorithm of Alder-Wainwright--. We shall analyze the interplay between the potential 
and the collisional dissipation typical of granular materials. In particular, we discuss the influence of the square-well 
interaction on the rate of energy dissipations in the same spirit of reference 7 -. It is well known, that in the homogeneous 
free cooling process, a system of inelastic hard spheres dissipates its kinetic energy at a rate proportional to the square 
root of the kinetic temperature, T, the so called granular temperature, and that T decreases in time as the inverse of 
square timei&. As we shall see, this picture is partially modified by the presence of a short range repulsive or attractive 
potential barrier. By treating collisions according to Enskog's equation, a generalization of Boltzmann equation to 
dense-fluid regime, we are able to make some predictions about the cooling behavior of the model. In addition we 
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consider its properties when it is kept in contact to a stochastic source of energy which balances the energy loss due 
to inelastic collisions. In this case, the system reaches a steady regime whose properties can be partly understood 
through a direct comparison with the properties of corresponding elastic system. 

The paper is organized as follows. Section II, describes the model we use and the main features of the simulations 
and technical details. Section III shows the thermodynamics of the elastic version of the model, in order to have a 
reference system to compare inelastic results. In section IV, an analytic estimate of granular temperature of the system 
is derived through a Boltzmann-Enskog approach. Section V illustrates simulation results of the inelastic model both 
in the cooling and driven regimes with a comparison to theoretical predictions. Finally section VI contains a brief 
discussion and conclusions. 



II. THE MODEL 



We consider N identical impenetrable rods of mass m = 1, size a = 1, positions Xi(t) and velocities Vi(t) constrained 
to move in a periodic domain of size L. They interact through a potential V(|xj — Xj\) consisting of a hard-rod part 
and a square-well potential, as shown in Fig. ^ Explicitly, we consider: 



V(x) = 



if x < a 
if a < x < ba 
if x > ba 



(1) 



where the parameter b defines the characteristic range of the interaction. The effect of a piece- wise constant potential 
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FIG. 1: Sketch of the interaction potential Eq.Q as a function of the interparticle distance x r , for 6 = 2 and a = 1. Solid line 
refers to attraction (e > 0) while dashed to repulsion (e < 0). 

V(x) amounts to a set of simple collision rules, similar to those involved in the dynamics of hard-rods. Several kinds 
of collisions between two neighbor particles may occur when their distance is: 

Axi(t) = Xi+i(t) - Xi(t) = ba 

If e > the following cases are possible: I) particles entering the well II) particles leaving the well, III) particles being 
trapped and rebounding at the inside square- well edge, because their relative kinetic energy is not sufficient to escape. 

On the other hand, if e < one has the cases: F) particles overcoming the repulsive barrier (m(Uj — w;+i) 2 > 4|e|) 
IF) particles descending the barrier, III') particles being repelled by the barrier, when m(vi — t>i+i) 2 < 4|e|. 

The post-collisional velocities in cases III and III' are given by: 

v'i = v i+ i 

W,Li = Vi. 



In the remaining cases the collision rule are found by requiring again the conservation of the total energy and total 
momentum at the edge of the square well-potential. If particles are entering (sj = sign(«i — i>;+i) > 0), while if they 
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are leaving the well (si < 0), and the collision rule reads: 



(Vi+v i+ i) (vi-Vi+i) 2 e 
Si\l : h Si — 

(2) 



(vi+Vi+i) (Vi~V i+ i) 2 6 
= Si\l h Si — 

4 m 



where pre-collisional and post-collisional velocities are indicated by unprimed and primed symbols, respectively. 
We finally consider the hard-core inelastic collision at 



which results in the transformation 



Axi(t) = x l+1 (t) - Xi(t) = a 
, l + a 

= V i+ l — (V i+ l - Vi) 

l + a ® 

v[ =Vi-\ — (Vi+l - Vi) 

where a indicates a constant coefficient of restitution and < a < 1 . 

Besides impulsive forces between particles, we shall also consider an external stochastic white noise force, whose 
role is to fluidize the system and balance the energy losses due to dissipative forces. The dynamics between two 
consecutive collisions is described by the following Langevin equation 

d 2 Xi(t) dxAt) . , 

where —m^dxi/dt is a viscous term and £i a Gaussian random force, with zero average and variance satisfying a 
fluctuation-dissipation relation: 

(&(*)&(*')> =2m 7 T<%<5(t-*') (5) 

with T proportional to the intensity of the stochastic drivingiii&. The damping term renders the system stationary 
even in the absence of collisional dissipation and physically can represent the friction between the particles and the 
container. Summarizing, the position Xi [i = 1, N) of the i— th particle evolves according to the equation: 

= -rnj^ + fc(t) + E /«(*) (6) 

j 

where fij indicates the resultant of impulsive forces between particles i and j. Since the dynamics of the model is 
mainly ruled by impulsive forces, Molecular Dynamics (MD) simulations make use of a collision driven algorithm. 



III. ELASTIC SYSTEM: EQUILIBRIUM PROPERTIES 

The elastic fluid, corresponding to the limit a — 1 in Eqs. @, serves conveniently as a reference system. Thus, we 
consider its equilibrium properties, that we shall compare to properties of the stationary inelastic system to build a 
theoretical approach valid in the region of moderate inelasticity. The equilibrium square-well fluid model is exactly 
solvable when the interaction range is restricted to first neighbors, i.e. b < 2, the exclude volume allows no more 
than two particles to experience the same potential well. In that case, the Gibbs free energy, G(P,T, N), can be 
derived using the isothermal-isobaric ensembl e) 1 ?' 21 ! 1 . Here, the partition function Y(P,T, N) is related to that of the 
one-dimensional canonical ensemble Z(T, L, N) by 



1 f°° 

Y(P T, N) = dLe~P pL Z(T, L, N) 



(7) 



where P is the thermodynamic pressure, A = h/y27rmkBT is the temperature-dependent de Broglie wavelength and 
Ao an arbitrary constant with dimension of a length. 

Following the existing literature^, the isothermal-isobaric partition function for N rods of length a can be written 

as 



Y(P, T, N) = A{_L [j'(e-PPo - e -/»«-) + e"^] }^\ (8) 



JV+l 
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The associated Gibbs potential 



G{P, T,N) = — ~ In Y(P, T, N) 



reads a part from a constant 

G(P,T,N) = (iV + l)|6aP+iln (/3PA) 

~ln[l + e^(e^ b - 1 ^ - 1)]}. (9) 

The equation of state, relating density, pressure and temperature is obtained by differentiating G with respect to P 
and defining the "volume" per particle, p^ 1 — L/N, of the system 

1 1 (b-l)a 

p = JP ~ 1 + ^Fi?p - l) • (10) 



that can be recast to the more familiar form 

(3P 1 



p 1 — poB 

with 

1 + c -/3PC-i)^( e -^-l) 



(11) 



Notice that 6=1 implies B — > 1, therefore the hard- rod pressure is straightforwardly recovered. 

In order to apply Enskog's kinetic approach, we need to compute the equilibrium pair correlation function in the 
thermodynamic limit. The pair correlation is defined as: 



PfJ 



( y ) = J2(S(x l+r - Xl -y)) (12) 



r=l 



To perform the average in eq. (|12|l . we represent the delta distribution (with I = 1 and q = r + I) by its Fourier 
transform 

f°° dk 

8{x q -x 1 -y) = J — exp[ik(x q - a?i - y)] (13) 
and write the average explicitly in terms of the Boltzmann weights f(x) — exp[— f3V(x)]: 



1 f 00 dk 



(5(x q -x x - y)) = z ^ - — J — e lky / dx q / dx q ^ • • • / dxi 







■f(L-x g )f(x g -x g _ 1 )---f(x 2 -x 1 )f(x 1 )e ik ^-^ . (14) 

Since the last expression has the form of an iterated convolution, one can obtain the desired average by means of 
standard Laplace transform method, a simple generalization of the method^ employed to compute Z(T, L, N) in the 
isothermal-isobaric ensemble. After a lengthy calculation, in the thermodynamic limit N — > oo and constant pressure, 
we get the series 

,, = f A§ (J3PYe-P p v 

P9\V) 2^1 (j- _ 1)1 [ e /3e( e -/3Po- _ e -/3P6<r) + e -pPbay- ^ 10 > 

where the coefficients can be written as 

^(v) = E(!"W-*('- + ( & - 1 )*)] 

k=0 ^ ' 

■[y - a{r + {b - l)fc)] r - X (l - e^f{e^) r ~ h (16) 
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where Q(x) is the unitary step functional. Notice that for a pure hard-rod system (e = 0, (3P = p/(l — per)), one finds 

oo r 

P9hr(y) = V t -Try- r^e(y - ra)(y - raf- 1 



■ exp 



1 — pa 



(y - rcr) 



(17) 



a result which agrees with the well known result by Zernike and Prins^S. From expressions ()15|l and , the value 
of the pair correlation at contact can be extracted explicitly: 



9&) 



(3P 



Since the thermodynamic pressure can be expressed in terms of g(x) through the virial equation 24 

— = 1 - [ dxV'(x)xg(x) 
P 2 J_ 0o 

after some rearrangements, the pressure reads 

^ = 1 + pcr[g(a) + bg(ba+) - bg(ba-)} 
P 

where the discontinuity of the potential at i = bo results in a jump of the binary correlation function: 

g(ba+) = g(a)e-^ b -^-^ 
g(ba-) = g(a)e-^ b -^ 



(18) 



(19) 



(20) 



(21) 
(22) 



with g(fe<T ± ) = lim,5_>o dibc ± <^)- These relations are consistent with Eq. (|l(Jf) and show that the pressure depends not 
only on the value of the pair correlation at contact, but also on its values at x = bo. 



IV. BOLTZMANN-ENSKOG EQUATION 

The system dynamics is determined by the combined effects of the heat-bath and interparticle collisions. Thus, 
the one-particle phase distribution function f(x,v,t), in the absence of external drift and large density fluctuations, 
evolves under the action of a Kramers operator associated with the interaction with the heat-bath 

^ 1 d 2 d 
K m(3 dv 2 ^ dv 

plus a collision operator, that we represent for the sake of simplicity as a Boltzmann-Enskog collision integral: 

df(v,t) 



dt 



C K f(v,t)+I(v,t). (23) 



Adapting to the present case arguments2^S& similar to those leading to the derivation of the SET (Standard Enskog 
Theory) we arrive at the following form of the Boltzmann-Enskog collision integral I(v,t): 

I(v, t) = I Q (v, t) + /+(«, t) + /_ (v, t) + I BS (v, t) (24) 

where the four contributions represent, respectively: 
• the inelastic hard-core collision 

I (v,t)=g(o) [ dv' [ dv"\v'-v"\ 



(25) 



■f(v')f(v"){S[v - l —^v' - ±±V] - 6[v - v"}} 
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the entering collision ((1/ — v") > 0) 

I+(v,t) = g(ba+) I dv' / dv"Q[{v' - v") 2 + 



4e, 



■5 



a — - + — 

4 m 



-5[v-v"]}\v'-v"\f(v')f(v") 



(26) 



• the escape collision ((v' — v") < 0) 

I-(v,t) = g{ba~) f dv' f dv"Q[(y' - v"f 

■U 



I] 

m 



(v' - v") 2 e 
4 m 



-S[v-v"]}\v'-v"\f(v')f(v") 



(27) 



• the elastic bound-state collision at Ax = ba 



I BS (v,t)=0 



which in one dimension vanishes and therefore can be omitted. 

We can apply the previous analysis to the theoretical description of the cooling process in the presence of the 
interparticle potential, under the hypothesis of spatial homogeneity. By integrating with respect to v the second term 
in the right hand side of Eq. I|25I27|I and approximating f(v) with Maxwellian distribution of temperature T g , we 
obtain the Enskog collision frequency at Aa; = a, and the frequencies of entering and escaping collisions: 



uo(po-) = (\v rel \)pg(a) 

u, + (pa) = {\v rel \)pg(ba+)[e(e) + 9(-e)e^] 

u-(pa) = (\vrei\)pg(ba-)[Q(-e) + e(e)e-^} 



(28) 



The expression for loq is formally identical to that obtained in the case of simple hard-rods without potential tail. 
It can be easily verified that the two factors containing the 0-functions are exactly the terms that compensate the 
asymmetry coming from expressions (|21() and l|22|l . Therefore, the rates become equal 



UJ + 



(\v rel \)pg{a)e-^ b -^ 
(\v re i\)pg(a)e-^ b -^-^ 



if e < 
if e > 



(29) 



and thus satisfy a detailed balance relation between entering and escape collisions. The presence of the potential is 
reflected in the modified value of the pair correlation at contact Eq. H18fl . By substituting (|iVez|) = 2( / 37rm) -1 / 2 we 
find the following expression for the collision time of the square well fluid: 



w 



P 



7Tml + e -PP(b-l)a^ 

For the sake of comparison the hard-rod Enskog frequency reads: 

J 



0e 



1) 



Uhr — 2\/ Phr- 



(30) 



(31) 



with (3P hr = p/(l — po). 

Interestingly, in an elastic system with repulsive forces, the ratio between the hard-core collision frequency and the 
entering/escape frequency reads 



J3P(f,-l)<T 



(32) 
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it suggests that, at high densities hard-core collisions dominate, because the pressure is a growing function of the 
density. Therefore, increasing the density amounts somehow to lowering the height of the effective potential barrier, 
i.e. the kinetic energy required to perform an elastic collision. 

We now consider, how the average kinetic energy of the inelastic system (a < 1) is dissipated. By multiplying 
Eqs. (|24() by v 2 and integrating over the velocity, we can compute the loss of kinetic energy due to collisions. Since 
only hard-core collisions dissipate energy we find that solely the process represented by Eq. I)25|l contributes to the 
evolution equation^ for the granular temperature T g (t) , 

^-(^"Vr (33) 

Notice that the expression (|30|l for wo, entering equation i|33|) . employs a pair correlation function g(a) extrapolated 
from its equilibrium value, where (3 has been identified with the inverse granular temperature, i.e. (3 — 1/T g . Moreover, 
the value of the pressure necessary to compute the frequency uj can be obtained numerically by inverting Eq. i|10|) 
for a given density of the system. The rate u>o decreases with increasing the repulsive barrier (e — * — oo) or when the 
temperature tends to zero. Consequently, as the system cools down, the dissipation rate will be much slower than the 
corresponding rate when e = 0. 



V. NUMERICAL RESULTS 



Whereas the equilibrium properties of the conservative system are analytically accessible, most of the properties 
of its dissipative version need MD simulations to be investigated. After resorting to numerical methods we shall 
compare their results with our theoretical estimates. At values of the restitution parameter a less than 1, the system 
is certainly not in thermodynamic equilibrium, but can achieve a stationary state when in contact with the heat-bath 
described by Eq. (JSJ). 

We shall consider the behavior of the system both in the cooling regime (7 = and T = 0) and in the stationary 
heated regime. The numerical methods employed have been briefly mentioned in the previous section and described 
in detail in papers 28-29 which the reader can refer to. 

In order to minimize surface effects and simulate an infinite system, periodic boundary conditions are imposed on the 
equations of motion. During each simulation run, we monitor the kinetic temperature, named granular temperature, 
T g , which by definition, is proportional to the average of the kinetic energy per particle: 



1 N 

-^mKt, 2 )-^) 2 ] (34) 



T — 
9 ~ N 

having chosen units in which fcs = 1. The values of pressure, instead, can be obtained from the virial formula 30 
properly modified for the present system: 

PL =l+ 1 + a m Z (35) 



NT a 2 t ob NT q 



where Z indicates the sum 



all coll 



t 00 is the observation time and ary = (xi — Xj) and = (vi — Vj) are, respectively, the separation and the relative 
velocity at the moment of collision. Both kinds of collisional events |scy| = u and \xij\ = bo determine an exchange 
of momenta among the particles. 



A. Cooling regime 



We consider, first, the properties of a system of N = 2000 particles evolving without the presence of heat-bath, 
thus no energy injection (T — 0) and no friction (7 = 0) occur. In the literature this situation is generally referred 
to as free cooling. The properties of this system with e = have been studied thoroughlyiSiiLi^il^ill and are well 
known. Due to the repeated inelastic collisions, the temperature T g decreases and after a short transient, lasting only 
few collisions per particle, T g displays the typical power law behavior t~ 2 , known as Haff's law. During such a regime, 
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the density remains homogeneous and the velocity distribution converges, from the initial Maxwellian, to a two hump 
function. As the system cools down, particles cluster into two "streams" at the outer edges the distribution and a 
bimodal velocity distribution emerges. 

Our MD simulations show that this scenario is modified by the the presence of potential tail (see Eq. QJ). Every 
MD run starts from an initial state characterized by N — 2000 particles with a Maxwellian velocity distribution of 
temperature T and uniformly distributed in space with no overlaps. During the dynamics, the grains spontaneously 
organize toward a state where the velocity distribution P(v) depends on the attractive or repulsive character of V(x). 
The behaviour of P(v) is clearly shown in figure O where later time distributions are characterized by a nearly 
Gaussian shape for e > and no longer Gaussian for e < 0. The attractive interaction has the effect to accelerate the 
dissipation, however the velocity distribution does not display the typical two-hump feature proper of the e = case, 
remaining a single peak function of shrinking width (see Fig[21left panel). 

On the contrary, when the potential is repulsive, only those pairs with velocities satisfying the condition (uj— Vi+i) 2 > 
4\e\/m may perform inelastic collisions. Such a selection mechanism is irrelevant when T g » e, i.e. the very early 
sta ges of the simula tions, however, it eventually leads to velocity distributions with small tails outside the region 
— y/\e\/m < v < y/\e\/m and almost flat inside (Fig. right panel) . Two small peaks can also be observed at 
v = ± y/\e\/m, likely, a reminiscence of the free-cooling two stream mechanism. 

Initial gaussian distribution 
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FIG. 2: Quenching of particle velocities observed at two stages of the cooling process of a system with attractive (e = 1, left) 
and repulsive (e = — 1, right) interparticle interaction. The histograms of the rescaled (dimensionless) velocity u = «ym/|e| 
are collected after n c hard-core collisions per particle have occurred. Simulations refer to a system of N = 2000 particles and 
parameters a — 0.99, pa — 0.002 and b = 2. The the Gaussian fits (dashed lines) are plotted for comparison. 



Equation (|30|l indicates that, under repulsive interaction, particles collide inelastically with an initial rate u>o oc yTg, 
that, as the system cools down, makes the crossover to the behavior ojq oc y/T^exp(e/T g ). Accordingly, fewer and 
fewer particle pairs will collide and the cooling slows down leading to a logarithmic decay in time of the temperature. 
However, this argument turns out to be incorrect. Indeed, Fig. [21 proves that the energy dissipation process occurs 
with a slower time decay than the prediction given by Eq. The direct comparison between theoretical and 

simulated dimensionless rate R = TT g /T is shown in the inset, where r is a proper time scale. The reason for such 
a discrepancy relies on the fact that the Maxwellian approximation for the velocity distributions, used to derive the 
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FIG. 3: Simulation results of the energy decay with time measured units r = ayrn/To, for e = (squares), e = — 1 (circles) 
and e = 1 (triangles) at effective density per = 0.002. Each point is the average over one-hundred trajectories of a system 
with TV = 2000 hard rods and initial temperature To = 10. Lines represent the analytical estimate from Eqs. (I33B and 13UL 
coherently with Eq. (1101 . Inset: Plot of the theoretical and numerical (dimensionless) dissipation rate R = tT versus the 
rescaled granular temperature (same symbols). 



rate expression (130(1 . fails as it seen from Fig. [2] (right) . With the actual shape of the distribution P(v), indeed, the 
system dissipates only a negligible fraction of the kinetic energy and undergoes an effective re-elasticization, implying 
that the non-Maxwellian character of P(v) is maintained up to the inelastic collapse. Our theoretical estimate of the 
collision frequencies works better at moderate densities, where dissipation can counterbalance the re-elasticization. 
As it suggested by Eq. I|32[l indeed, The pressure exerted by the dense surrounding fluid on two colliding partners 
may overwhelm their repulsion, so that they will experience frequent hard-core collisions, i.e. loq » omega±. 



B. Driven regime 



The scenario changes when the system is coupled to an heat-bath at temperature T. A steady regime, characterized 
by almost constant granular temperature and pressure, is attained. As already done for the cooling regime, we can 
derive an implicit relation for Tg&2& by multiplying both sides of Eq. I|23[) by v 2 and integrating with respect to v, 

T g =T (37) 

The variation of T g with density, given by the numerical solution of Eq. I|37|l ■ is compared in figure 0] with the results 
of MD simulations. The agreement between theory and numerical experiments is satisfactory for the three possible 
cases: attractive, repulsive and vanishing inter-particle interaction. The virial formula (|35fl is employed to compute 
the pressure of the system by averaging over different MD runs. The simulated pressure values are plotted, in figure 
03 together with those obtained by a self-consistent solution of formula (|20|l with the appropriate replacement of the 
heat-bath temperature T by the granular temperature T gi Eq. (|37fl . 

The use of formula (|20l) implicitly assumes that the pair correlation function for the inelastic model maintains the 
same functional dependence as its equilibrium counterpart. Such an hypothesis can be checked by measuring during 
MD runs the three collision frequencies luq, lo+ and uj— and comparing them with their theoretical prediction. The 
behaviour of these quantities with the dimensionless variable pa is reported in Fig. for both attractive and repulsive 
interactions. Even in the inelastic case, one observes that the ratio, u>q/uj±, between frequencies of dissipative collisions 
and barrier crossing increases with density from the value 1, observed in a very diluted system, as shown in Fig. 
This is very consitent with the prescription provided by formula (|32fl that, hard-core collisions, in this model, become 
dominant events at higher densities. 

In the case of barriers (Fig. Eji), the theoretical frequencies agree fairly with those extracted from simulations. 
However, some discrepancies arise when particles may mutually attract (FigJHJa), even though the overall trend of the 
frequencies with the particle density is correctly captured by the theoretical predictions. The differences induced by 
inelasticity become more evident in Fig. El where we plot the theoretical and numerical pair correlation functions g(y) 
at the value pa = 0.05 for repulsive c) and attractive d) particles. Again, the theory fits faithfully the simulations for 
the system with repulsive interactions, while, for attracting particles, the simulated g(y) deviates of about a factor 
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FIG. 4: Dependence of the ratio between granular T g and bath temperature T = 10 on the density, for a system with repulsive 
e = — 10 (squares), vanishing e = (circles) and attractive e = 10 (triangles) interparticle interaction. Points indicate the 
average over a set of 10 4 samplings in a single MD run of duration t ma x = 10 4 . Solid lines refer to theoretical temperatures 
extracted from the numerical solution of Eq. 13711 . The number of particles is N = 2000, the remaining parameters are chosen 
as a = 0.9, 7 = 0.2 and 6 = 2. 
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FIG. 5: Pressure of an inelastic system, rescaled to the equivalent ideal gas pressure {Pid ~ pT), as a function of pa in the case 
of e = — 10 (squares), e = (circles) and e = 10 (triangles), for a system with the same parameters as in Fig. gl Solid lines are 
the corresponding analytical values obtained according to formula 1201 1. 



two from its estimate in the interval a < y < ba. Figures (FigEK) and (FigJHt>), on the contrary, indicates clearly 
that our theoretical approach perfectly describes the functions g(y) of the elasitc system with both attractive and 
repulsive forces. 

The difficulty encountered by the theory to fit some regimes of the system with attracting inelastic particles can 
be ascribed to the different effects that repulsive and attractive interactions induce on the inelastic system. The 
former basically entails a system re-elasticization which may favor homogeneous particle distributions, while the 
latter enhances the frequency of inelastic collisions leading to clustering. The relevant physical parameter controlling 
the system behaviour is the ratio |e|/T fl , thus a good evaluation of T g is crucial to make accurate theoretical predictions. 
Our estimate of the granular temperature T g , from Eq. JHZJl, is based on the assumptions of spatial homogeneity. For 
repulsive interactions (barriers) the homogeneous state occurs, while for cohesive interactions, particles, under specific 
conditions, can easily cluster making the system inhomogeneous. If this happens, the single observable, T gi does not 
describe properly the kinetic state of the system and, in addition, its estimate from Eq. (|37l) is incorrect since that 
formula neglects local temperature fluctuations. 

The deviations of the theory from simulations become less pronounced as 7 increases, and the reliability of the 
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FIG. 6: Collision frequencies at particle separation xtj = a and Xij — bo as a function of pa. Figure (a) shows the repulsive 
case (e = —10) while figure (b) refers to attractive interaction (e = 10). Lines correspond to the theory from Eqs. The 
remaining parameters are as in Fig. 2] 

theoretical approach can be quantified by the integrated difference between numerical, g n (y), and theoretical, gt(y), 

A<? = -/ dy[g n (y) - g t (y)] (38) 

for various values of 7, but e/T — const. The dependence of Ag on 7 reflects the fact that the response of the fluid 
to the action of the heat-bath is faster as 7 — > 00 and thus erases more rapidly the memory of inelastic collisions. 
Within this limit one recovers the behavior of the elastic system. 

VI. CONCLUSIONS 

In this paper we have investigated both theoretically and numerically the influence of a finite range interparticle 
interaction on the behavior of a one-dimensional inelastic hard-rod system. Forces and interactions, whose range 
is larger than the size associated with the excluded volume constraint are often present in many realistic granular 
materials. In the specific case, we have chosen a square well potential to model attraction and a square barrier to 
model repulsion. These simple shapes, in the case of undriven system, still enable a computer implementation of 
the particle evolution in terms of a collision driven molecular dynamics. In fact, simple transformations describe the 
instantaneous changes of velocities when the separation between two particle corresponds to the two characteristic 
ranges of the potential. We first analyzed how the interplay between these finite range forces and inelasticity modifies 
the cooling scenario with respect to the free inelastic system. We found that in the case of repulsive barriers the 
temperature decay becomes slower than Half's 1/i 2 power law and eventually reaches a regime where the system is 
nearly elastic. In the case of attractive wells, instead, the granular temperature is lost faster than an inverse time 
power law. 

Secondly, we studied the behavior of the stationary regime obtained through a stochastic forcing of the system. 
The steady state has been analyzed via MD simulations and theoretical approaches based on the direct comparison 





12 



8.0 



6.0 



r 2 - 4.0 



2.0 



CO 

o 

CO 



o.o 



T g =0.Mel 

T =0.5lel 

g 

T =lel 




0.1 



0.2 

pa 



0.3 



0.4 



FIG. 7: Collision ratio (Eg. 1321 . as a function of pa at temperatures T g — \e\, T g = 5|e| and T g = 10|e|, for a driven system of 
N = 2000 particles with inelasticity a = 0.99 and a barrier of sizes b — 2, e = —1 in contact with a bath of viscosity 7 = 0.2. 
Each point is the average over a single run of duration t max = 10 4 , sampled every t Bamp i e = 1.0 time units = I/7. 

with the elastic counterpart of the system whose equilibrium properties are well understood. Our results show that, 
in the dense limit, particle spatial correlations are relevant and modify the collision rate, the excluded volume of the 
other particles enhances the probability that two particles are at contact and thus reduce the repulsive barrier. The 
theoretical approach we have attempted remains a reliable approximation for the behaviour of the dissipative system 
at not too small densities while it is correct for the elastic system at every density. 
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FIG. 9: Parametric plot of the cumulated difference Ag between theoretical and simulated correlations function (see Eq. 
versus the friction coefficient 7, for an inelastic system with N = 2000 rods and parameters a — 0.9, T = 10, e = 10, pa = 0.05 
and b — 2. 



